if (patchDEMID > -1)
{
    fixedPotentialIDList.clear();
    fixedPotentialValueList.clear();
    volScalarField cellFlux(fvc::div(phi));
    forAll(topCellID,celli)
    {
        label currentCell = topCellID[celli];

        if(h[currentCell] >= distanceToDEM[celli])
        {
            if (cellFlux[currentCell] < 0)
            {
                fixedPotentialIDList.append(currentCell);
                fixedPotentialValueList.append(distanceToDEM[celli]);
            }
        }
    }

    // Display number of seepage only for serial run (wrong value for parallel
    label nSeepageCells = fixedPotentialIDList.size();
    reduce(nSeepageCells, sumOp<label>());
    Info << "Number of seepage cells = " << nSeepageCells << endl;
}
